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We investigate the possible use of Bose-Einstein condensates of diatomic molecules to measure 
nuclear spin-dependent parity violation effects, outlining a detection method based on the internal 
Josephson effect between molecular states of opposite parity. When applied to molecular conden- 
sates, the fine experimental control achieved in atomic bosonic Josephson junctions could provide 
data on anapole moments and neutral weak couplings. 
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Since the discovery of the electroweak unification, par- 
ity violation (PV) has been one of the outstanding prob- 
lems in atomic and molecular physics. Although the suc- 
cess of the standard model (SM) of elementary parti- 
cles is extraordinary, the search for new physics is still 
carried out at high- and medium-energy particle collid- 
ers. On the other hand, very sensitive low-energy probes 
can also be implemented to test SM predictions As 
compared to large accelerators, atomic experiments offer 
a complementary approach to study PV. In particular, 
nuclear spin-dependent PV effects focus on two causes. 
One is the electroweak neutral coupling between elec- 
tron vector- and nucleon axial- vector currents. This can 
be parametrized by two constants, C2„/d, describing the 
electron couplings to up/down quarks. These constants 
are the most poorly characterized parameters in the SM 
[1]. Thus, precise measurements of C2u/d are potentially 
sensitive to new physics beyond the SM Q. Another 
cause of spin-dependent PV is the weak nucleon-nucleon 
interaction which gives place to the anapole moment 
0, [Ij. This parity violating magnetic moment, which re- 
sults from the chirality acquired by the nucleon current 
0, couples to the spin of the electron. As claimed by 
Haxton and Wieman 0, the most practical strategy for 
studying the effects of Z'^ exchange between hadrons is 
the investigation of the parity- violating nucleon-nucleon 
interaction. One of its signatures is the resulting anapole 
moment, which has only been measured in ^^"^Cs Q. 
These high precision measurements allow for a better 
determination of the meson-nucleon parity violating in- 
teraction constants, with standard reference values given 
by the Desplanques, Donohue and Holstein (DDH) best 
values 0]. As there are large uncertainties in these pos- 
sible values (the DDH reasonable ranges), it is therefore 
important to determine the weak coupling constants ex- 
perimentally. 

These and other important tests of the SM could be 
performed by using the ultra-high resolution which is po- 
tentially available in cold molecule experiments. In par- 
ticular, it is known that the internal structure and spe- 
cific properties of diatomic molecules can enhance the 
violation of some discrete symmetries, as compared to 
their atom counterparts [lo| . In addition, heteronuclear 
molecules exhibit dipole-dipole interactions that make 
them attractive candidates for use in quantum simula- 
tions of condensed-matter systems and in quantum 



computation [12| . Thus, experiments involving ultracold 
heteronuclear diatomic molecules will be valuable to a 
large portion of the physics community. The fact that re- 
cently both homonuclear [l^l and heteronuclear [T5| 
ultracold molecules have been produced in the rovibra- 
tional ground state, suggests that the type of experiments 
envisaged here may not lie too far in the future. 

The total nuclear spin-dependent parity violation effect 
is given by the Hamiltonian [l| 



Hpv = —Kyj-j—^o. ■ Ip(r), 
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where Gp is the Fermi weak constant, I is the nuclear 
spin, K = {I + l/2)(-I)^+i/2-'^ / is the orbital angu- 
lar momentum of the unpaired nucleon, a is the vector 
whose components are the Dirac matrices acting on the 
electron spinor, and p(r) is the nuclear density. The re- 
term is 
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where Ka stands for the anapolar term, K2 is due to 
the electroweak neutral coupling between electron vector- 
and nucleon axial-vector currents and kq represents the 
interaction due the nuclear weak charge perturbed by 
the hyperfine interaction [T|. We note that the nuclear 
anapole moment is the dominant source of PV in atoms 
with atomic number ^ > 20 since it increases as 
Ka ^ g J^l'^ ; where g is the strength of the weak interac- 
tion between the unpaired nucleon and the nuclear core 
(the Kq term also increases as A^l'^ but the numerical 
coefficient is very small). Let us recall that there have 
been proposals for measuring the anapole moment in 
the ground state of heavy alkali-metal atoms [l^ . Elec- 
troweak calculations show that, for a heavy nucleus such 
as, e.g., Rb, Kq/kq = 20 and Ka/n^ = 5 |17l |. 

Concerning molecules, we note that the levels of op- 
posite parity are four to five orders of magnitude closer 
in diatomic molecules than in atoms, which results in 
stronger mixing of those levels by PV The fact 

that the near-degeneracy of a pair of opposite parity 
levels enhances the level mixing caused by the weak in- 
teraction can be exploited for precision measurements, 
as shown by the Berkeley group with atomic dyspro- 
sium [131 ■ We note that this enhancement mechanism 
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does not work for homonuclear diatomic molecules [19|. 
Moreover, this energy interval between opposite parity 
levels can be made even smaller by applying an external 
magnetic field. In addition, and contrary to the case of 
anomalously close levels of opposite parity in rare-earth- 
metal atoms, the natural linewidths of all the levels in 
the electronic ground state of a molecule arc negligible 
[20| . Most often, the ground state of a molecule which 
has an odd number of electrons is either the or the 
^Ili/2 state. However, since the rotational structure of 
the lower energy levels is simpler in the former case, we 
will concentrate on diatomics, which are the molecu- 
lar equivalent of alkali-metal atoms. 

Almost all the molecules that have as their ground 
state pertain to the Hund's b coupling scheme. Follow- 
ing Ref. poj . L (electronic orbital angular momentum) 
is coupled to R (rotational angular momentum of the nu- 
clei) to form N, and N is coupled to S (electronic spin) 
to form J. Their corresponding projections on the in- 
ternuclear axis, n, are, respectively, m/,, mjv, ms and 
mj. Including nuclear hyperfine interactions by the cou- 
pling J -|- I = F, one obtains the spin-rotation-hyperfinc 
Hamiltonian: 

i/srh = SN2+7N-S + 6I-S + c(I-n)(S-n), (3) 

where 7 is the spin-rotation constant and b, c are hyper- 
fine constants. 

In most cases of interest, B ^ 7, &, c, so is a 
good quantum number with eigenstates of energy En ~ 
B N{N + 1) and parity P = (—1)^. Since one can use a 
magnetic field to Zeeman shift consecutive states of op- 
posite parity {N = 0, 1), and the magnetic field necessary 
to overcome the rotational energy is large enough (about 
1 T when the rotational constant is i? ~ 1 GHz), one can 
work in the decoupled basis \N,'m]\f)\S,ms)\I,rni) . 

Within the subspace of rotational-hyperfine levels, the 
parity violating Hamiltonian ([1]) can be written effec- 
tively [13 

ffpv ^WpK{nxS)-I, (4) 
where Wp can be calculated using semiempirical methods 

The opposite-parity levels 

\A) = |0,0)|l/2,m5)|/,TO/) 

|B) = |l,m'^)|l/2,-ms)|/,m;) (5) 

can be mixed by iJpv when the quantum number mp = 
niN + ms + mj is identical for both states lA) and \B) 
p^ . Thus, once Wp has been calculated [2l| together 
with the matrix elements of the operator (n x S) • I be- 
tween opposite-parity states, we seek to determine k by 
measuring 

iW = KWp{A\{nxS) -IIB). (6) 

Importantly, time reversal invariance ensures that iW is 
purely imaginary. We note that this scheme is currently 



being designed to measure iW using a Stark-interference 
technique |23| . 

Here we focus on the effect of iJpv in an internal 
bosonic Josephson junction (BJJ) between two rotational 
states in a gas of Bose condensed diatomic molecules, 
achievable with foreseeable advances in molecular laser 
cooling. In this case, PV effects can be manifested in 
the population difference and in the relative phase be- 
tween the two modes, which can be directly measured 
both in the Rabi and Josephson regimes [IJ]. The 
description of the Gross-Pitaevskii (GP) dynamics re- 
duces, at low energies, to a nonlinear two-mode equation 
for the time-dependent amplitudes, ipi{t) = Ni{t) e*^' 
{i = A,B), where Ni and 9i are, respectively, the con- 
densate molecule number and phase in state i. When the 
two condensates are coherently linked, we write the wave 
function as 

|*(t))=VA(i)l^)+^B(0|S>, (7) 

where iVA(i) + ^s(i) = Nt- 

Including i?pv, the amplitudes and phases obey the GP 
equation 
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where Ei is the effective energy (which includes the mean- 
field contribution), and J7 is the highly tunable Rabi fre- 
quency connecting the two modes |25l - [27| . 

Defining the fractional population imbalance, z{t) = 
[NA{t) - NB{t)]/NT, the relative phase, (/>(*) = 6 Ait) - 
Osit), and rescaling to a dimensionless time J7t — ?> t, we 
obtain 

i = — \/l — z'^ (sin (f> + w cos 0) 

<j) = Az + e (9) 
2; 

H -^=^= (cos (f> — w sin </>) , 

\/l — z'^ 

where e is the effective detuning, A = UNt/^, w = 
2W/il and U encompasses the various interactions. We 
note that these equations represent a generalization of 
the well known pendulum equations that describe the 
classical dynamics of two weakly coupled Bose-Einstein 
condensates [l^H^ . They clearly show the PV effect in a 
molecular BJJ, with the parity violating contribution W 
competing with the usual coupling matrix element il/2. 
We conclude that electroweak effects will have measur- 
able consequences on the macroscopic dynamics of the 
condensate. 

Equations ^ derive from the Hamiltonian 

Az^ 

H = — — h — vl — (cos (j) ~ w sin 0) , (10) 

which describes a parity-violating non-rigid pendulum, 
where the parity operation in this representation changes 
(j) — > —(j). Within this mechanical analog, PV is re- 
fiectcd in the shifted equilibrium position of the pendu- 
lum, which now is — arctaniy due to the s'mcj) term 
in (Uni). 
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For e = and negligible interactions (A 
oscillations have the modified Rabi frequency 



0), the 



n = 0^1 + (11) 

In this simple case, the population imbalance evolves as 



z{t) = 



1 + 



sm{nt + A2), 



(12) 



where Ai = 1 



zoVi^ + w^)/Ai 



with zq = 



and A2 = 
;(0). The evo- 



sm 

lution of (j){t) is also very sensitive to iW but cannot be 
written in closed form. 

Deviation from the standard non-rigid pendulum dy- 
namics is shown in the time dynamics of both the popu- 
lation imbalance and phase difference (Fig. [T]). 





Figure 1. Time evolution of the population imbalance and 
phase difference for A = and w = (solid), 0.3 (dotted) and 
1.0 (dashed). We have taken zq — 0.6 and (j>o — n. 



an identical response from all molecules, which should 
improve the accuracy of the experiment. Following Ref. 
j2g | , we estimate that the condensate response to a spa- 
tially varying Rabi frequency in the z axis is rigid when 
El ~ Eq > |r2'|(l|z|2), where V,' stands for the gradi- 
ent of the Rabi frequency, and 1 and 2 refer to the two 
lowest states of the trap. Taking typical values for trap- 
ping frequencies of 27rx5 Hz and trap sizes of about 25 
fini, we obtain that the condensate response is rigid when 
<2ttx 0.1 Hz /im-i. 
In the strongly nonlinear regime, macroscopic quan- 
tum self trapping (MQST) is also present [2^. In fact, 
the value zq = is inaccesible at any time if 



A> Ac 



A, 



w 



1 — Zq sin ( 



(15) 



where c/jq = </>(0) and Ac = ^ (^■\/T~~z^ cos 0o + 1^ is the 

critical parameter for MQST in absence of PV effects. In 
Fig. [5] we show, for A = Ac = 1.11, the comparison of 
parity-even and parity-odd effects in the Fourier trans- 
form of the population imbalance, noting that the PV 
effects are easily distinguishable from their parity con- 
serving counterpart. We notice that the harmonic con- 
tent of z{t) is quite sensitive to iW for large A. In the 
inset we explicitly show the effect of PV in the time evo- 
lution of the population imbalance. We note that, due 
to the difference between Ac and Ac, the onset of MQST 
depends on the value of w. 
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In the presence of interactions [28| , the harmonic limit 
(|z|,|0 — ^ 1) displays sinusoidal oscillations of z{t) 
with a modified Josephson-Rabi frequency 



wjR, = O (1 + A) 



1/2 



(13) 



where A = A/VT + i/^ is the PV-renormalized interac- 
tion parameter. The linearized Eqs. ^ read 







= 0. 



(14) 



We conclude that in the harmonic limit, PV effects 
change the Josephson-Rabi frequency. In the non- 
interacting limit, the physics is formally identical to that 
of a non-degenerate gas of molecules, if the external fields 
are uniform. In the presence of inhomogeneities, an im- 
portant advantage of the condensate is that it permits 



Figure 2. Time evolution of the population imbalance (inset) 
and its Fourier transform for A = Ac = 1.11 and w = (solid) 
and 0.2 (dotted). 

Once an appropriate control over the internal states 
of diatomic molecules is achieved, one also needs a pre- 
cise measurement of the conjugate variables [24| to 
measure PV effects in a molecular BJJ setup. Hereafter 
we focus on the molecule ^^Rb'^Yb (with even A). Yt- 
terbium is so far one of the few elements besides the 
alkali metals that can be trapped and condensed in opti- 
cal traps [lOl ■ In addition, as Yb is a closed s-shell atom 
with zero nuclear spin, the only active electron in the 
molecule is the Rb s-electron. Thus, from the point of 
view of the electroweak interaction, the Yb atom plays a 
passive role. 
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First we estimate the value of k. ^^Rb has 7 = 3/2 
with an unpaired proton in the p state. Thus, k 



100 ms. Following Ref. [3J we take / = 20 and get that 



= the effect of Hp^ could be experimentally detected within 



+ |k2 — ^'^Q- To get an insight on the order of mag- 
nitude expected for the parity violating effect, we con- 
sider K ~ 0, [l3|- Taking gp — 5 we get IOOkq — 
27. Given the proportionality between matrix elements of 
different weak interactions, we estimate the value of Wp 
using previous molecular electronic structure calculations 
of the electron electric dipole moment of RbYb [3l[ , get- 
ting Wp ~ 36 Hz (for comparison, we note that RaF has 
one of the largest values for Wp predicted so far, which is 
in the range of kHz [s^l). Thus, we conclude that the ex- 
tra coupling that has to be accounted for in a RbYb-BJJ 
is, including only the anapolar term, ~ 4Hz. If we 
include also the contribution of K2 — 0.05 and kq ~ 0.01 
p7| . we obtain W ~ 5 Hz. 

The Rabi frequency $7 is a highly tunable parameter 
because of its linear dependence on the applied dc electric 
field. We note that RbYb has a typical value for the 
electric dipole matrix element of 1 kHz cm V~^. Thus, 
for instance, a value O ^ 100 Hz can be obtained using an 
electric field ~ 0.1 V cm^^ . Due to the small value of this 
matrix element, dipole-dipole interactions are negligible 
compared to both and W. In addition, we note that 
Zecman degeneracy can be controlled up to 10"^ Hz for 
RbYb due to the high level of magnetic field control (few 
tens of /iG) recently achieved in a cold atom experiment 



The curves z{t) [Eq. ([T2|) ] and (f>{t) (which lacks a 
closed form but can be calculated numerically with arbi- 
trary accuracy) depend onw = 2W/fl through their fre- 
quency and amplitude. If we focus on the Rabi regime, 
where interactions are negligible, and assume that we es- 
timate the value of w only from the measurement of the 
oscillation frequency, the resolution obtained for IF ~ 5 
Hz is known to be AW ^ f / {t^/Nt), where r is the 
coherence time and / is numerical factor taking into ac- 
count recent experimental capabilities in the absorption 
imaging of ultracold molecules. Wc take Nt — 4 • 10** as 
the total number of condensed molecules and assume r ~ 



a relative precision of ~ 20%. We point out that this 
improvement in the accuracy of AVF is due to the long 
coherence times of trapped cold molecules, as compared 
to those of a molecular beam experiment j32| . Concern- 
ing this coherence time, an important point is the role 
of losses. For RbYb wc expect a loss time of about 100 
ms [1^. This time should be larger than both the Rabi 
and PV times, and W~^, for this proposal to be vi- 
able. The value of ft is not a problem since it is highly 
tunable (or even zero). For RbYb, taking ~ 5 Hz, 
this condition can hardly be satisfied. However, we note 
that current estimates of Wp are still rather crude and 
the real situation might turn out to be more favourable. 
On the other hand, molecules such RaF, HgH or HgF are 
known to have much larger values of Wp |32| . 

Understanding the final sensitivity that could be 
reached by using a molecular B JJ requires a careful study 
of systematic effects. As noted in Refs. [il,!!^' care has 
to be taken of stray electric fields that could mimic the 
PV effect. Due to the small value of the electric dipole 
matrix elements, it suffices to control stray fields within 
a feasible accuracy of 0.1 V cm~^. Thermal effects can 
be neglected if the temperature is smaller than a typical 
rotational level splitting (~ mK). However, temperature 
can be bigger than W and fi, thanks to the fine control 
of the initial state. 

We finally note that in the proposed setup a large num- 
ber of curves depend on a few parameters. This intrinsic 
redundancy will permit not only a consistency check of 
the underlying physics but a precise measurement of the 
anapolc moment in a given molecule. The eventual mea- 
surement of this moment in a variety of molecules may 
ultimately reveal the strength of the poorly known elec- 
tron couplings to up and down quarks. 
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